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Abstract:  The  adhesion  interaction  of  an  Ir  tip  with  Ir  and  A1  surfaces  has  been 
investigated  using  the  scanning  tunneling  microscope.  QuaOtative  diffe.rences  in  the 
adhesion  characteristics  have  been  observed.  A  duster  model  based  on  effective 
pair-potentials  has  been  set  up  to  understand  adhesion  in  atomic  scale  contacts. 
The  model  calculation  provides  a  semi-quantitative  explanation  of  the  experimental 
data.  In  particular  it  is  shown  that  subtle  differences  of  the  atomic  interaction 
potentials  have  a  significant  effect  on  the  contact  mechanics. 
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The  study  of  metallic  adhesion  has  a  long  history  yet  our  understanding  of 
atomic  scale  processes  is  still  sketchy  [I].  Valuable  insight  has  been  provided  by 
recent  molecular  dynamics  simulations  [2]  but  very  little  is  known  experimentally 
about  the  short-range  chemical  adhesion  interaction  which  Ls  responsible  for  the 
metallic  bonding  in  the  contact  zone. 

In  this  letter  we  focus  on  the  interaction  between  a  very  sharp  metal  tip  and 

a  flat  metal  surface.  The  interaction  is  investigated  experimentally  using  an 

ordinary  scanning  tunneling  microscope  (STM)  with  a  cantilever  spring  (CB)  as 
$ 

sample  stage.  The  tip  of  the  STM  serves  simultaneously  as  a  force  probe  and  for 
tunneling.  Tne  latter  is  exploited  to  establish  a  well-defined  distance  between 
probe  tip  and  sample.  Direct  measurement  of  the  tip-sample  force  is  difftcult 
because  stiff  CBs  must  be  employed  to  prevent  instabilities  of  the  gap.  The  force 
gradient)  on  the  other  hand,  is  readily  detected  by  means  of  a  djmamic 
measurement  [31.  The  basic  principle  is  that  the  force  gradient  couples  the  tip  to 
the  CB  causing  the  resonance  frequency  of  the  CB  to  shift  in  proportion  to  VF. 

Experiments  were  conducted  under  ultrahigh  vacuum  conditions 
{p  <  3  X  10~^^  mbar).  Tunneling  tips  were  made  of  mechanically  sharpened  Ir 
wire.  Final  tip  forming  was  performed  in  siiu.  Ir  samples  were  prepared  by 
extensively  sputter  cleaning  the  surface  of  CBs  made  of  po[ycr5'’stalIine  Ir  sheet 
metal  Typical  length,  width,  and  thicknes.s  of  the  ^ams  were  5  mm,  0.5  mm  and 
50  fxm,  respectively,  yielding  a  resonance  frequency  of  =il.5kHz.  The  spring 
constant  —  typically  on  the  order  of  100  Nm~^  -  was  adjusted  by  positioning  the 
tip  between  the  clamped  and  the  free  end  of  the  CB  which  had  a  nominal  spring 
constant  of  =^50  Nm“^  AI  samples  were  prepared  by  condensing  thin  films 
(sslOOO  A)  onto  a  substrate  Ir-CB  at  room  temperature  and  a  base  pressure  of 
/?  <  5  X  10“^°  mbar  [4]. 
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The  interaction  force  gradient  was  measured  as  a  function  of  tip-sample 
distance,  z,  defined  as  the  displacement  of  the  tip  with  respect  to  a  reference 
position  characterized  by  a  specific  value  of  the  tunnel  resistance, 
Ry  —  10‘ Q  /  10^  Q.  for  the  Ir  and  A1  samples,  respectively.  With  the  tip  at  the 
reference  position,  the  feedback.  loop  which  adjusts  for  a  constant  tunneling 
current  was  interrupted  and  the  tip  was  ramped  towards  the  surface  at  a  rate  of 
The  tunneling  current  and  the  resonance  frequency  were  simultaneously 
recorded  during  approach-  In  order  to  prevent  accidental  touching  of  the  surface 
the  tip  was  rapidly  retracted  to  the  reference  position  as  soon  as  the  tunneling 
resistance  dropped  below  ~I00kn.  After  each  approach  cycle  the  feedback  loop 
was  turned  on  for  100  ms  to  compensate  for  possible  drifts. 

The  results  of  the  experiments  are  summarized  in  Fig.  1.  Each  enrve 
represents  an  average  of  64  approach  cyclas.  Interaction  force  gradients  measured 
on  the  Ir  sample  are  negative  and  rapidly  increase  in  magnitude  as  z  approaches 
—2.6  A  which  is  the  maximum  tip  excursion  allowed  by  the  tunnel  resistance 
constraint  The  observed  adhesion  characteristics  are  roughly  consistent  with  the 
concepts  of  universal  adhesion  developed  by  Rose  et  al.  [5,6].  In  a  separate 
experiment  the  tunnel  resistance  limit  was  ignored  in  order  to  probe  the  contact 
regime.  At  zoi  —2.7  A  the  tunnel  resistance  abruptly  drops  by  one  order  of 
magnitude  and  subsequently  levels  off  at  c^lO  kii.  Wc  were  not  able  to  measure 
interaction  force  gradients  reliably  in  this  regime.  In  addition  hysteresis  effects 
were  observed  which  point  to  irreversible  changes  of  the  atomic  structure  of  the 
tip  or  the  sample  surface.  The  tunnel  resistance  curve  for  r<— 2.6A  is  therefore 
complemented  by  one  representative  example.  Adhesion  as  well  as  tunnel 
rcsi^nce  cnaracteristics  for  the  Al  sample  are  distinctively  different.  At  large 
tip-saraple  distance  the  interaction  force  gradient  is  negative,  similar  to  the  Ir 


sample,  but  its  magnitude  docs  not  increase  monotonically  with  decreasing  gap 
width.  Instead,  a  minimum  is  observed  and  eventually  the  force  gradient  even 
becomes  positive.  Exponential  tunnel  resistance  characteristics  are  observed  in  the 
regime  of  negative  interaction  force  gradients.  Simultaneously  with  the  force 
gradient  becoming  positive  the  tunnel  resistance  curve  levels  off  and  approaches  a 
value  of  ~I00kI2.  No  indication  of  a  discontinuity  in  the  the  tunnel  resistance  at 
small  gap  widths  can  be  discerned.  However,  the  tunneii.ng  current  was 
significantly  less  stable  for  the  AI  samples  than  for  the  Ir  samples. 

To  understand  adhesion  in  atomic  scale  metal  contacts,  we  have  set  up  a 
cluster  model  to  calculate  the  forces  between  a  metal  tip  and  a  flat  metal  surface. 
The  surface  is  that  of  a  three-layer  duster  of  hexagonal  symmetry  with  a  total  of 
21 1  atoms.  The  tip  consists  of  a  base  layer  of  4S  atoms,  a  second  layer  of  12 
atoms  and  a  single  additional  atom  on  top  of  that.  Clusters  of  this  size  cannot  be 
treated  ab  initio  at  present,  so  the  setni-empirical  embedded  atom  method  is  an 
attractive  alternative  [7].  Unfortunately,  the  relevant  parameters  for  Ir  have  not 
been  determined  as  yet  We  therefore  resort  to .  a  model  in  which  the  total  energy 
is  calculated  by  summing  nearest  neighbor  two-body  interactions.  For  this  we 
could  use  the  interaction  between  two  isolated  atoms,  calculated  by  some 
first-principles  method.  Summing  such  (pure)  two-body  interactions  to  obtain  the 
energy  of  large  clusters  would,  however,  ignore  all  many-body  effects  in  metals. 
We  therefore  dctcrmuie  effective  nearest  neighbor  interactions  by  calculating 
smaller  clusters  consisting  of  seven  atoms  arranged  in  two  parallel  triangles  with 
the  seventh  atom  in  between.  To  get  the  effective  interaction  energy  between  two 
lilte  atoms,  we  assume  all  distances,  d,  in  the  cluster  to  be  the  same  and  define 


(I) 
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where  is  the  total  electronic  energy  of  a  cluster  consisting  of  n  atoms.  To  get 
the  effective  two-body  interaction  between  two  different  metals,  we  take  a 
tetrahedral  cluster  of  the  first  metal  interacting  with  a  triangular  cluster  of  the 
second  metal  at  a  distance  d  from  the  tip  of  the  tetrahedron  and  define 

=  y(£7  -  4''''’  -  4"”)  ■  (2) 

The  bond  lengths  in  the  individual  clusters  are  adjusted  to  minimize  the 
respective  energies,  thus  obtaining  their  respective  equilibrium  geometries.  For  the 
calculation  of  (1)  and  (2)  we  have  used  the  ASED-MO  method,  a  semi-empiricai 
tight-binding  model  with  matriT:  elements  calculated  in  an  extended  Huckel  scheme 
[8].  In  Fig.  2  we  plot  (1)  and  (2)  rescaling  the  energy,  V,  and  the  distance,  by 
the  respective  values  at  the  potential  minimum.  The  Ir-Tr  potential  (solid  line)  falls 
off  exponentially  with  a  decay  constant  as  expected  from  universality  theory  [5]. 
It  is  important  for  our  present  study  that  the  Al-AI  potential  (dotted  line)  is 
different.  In  particular,  it  falls  off  much  more  slowly  with  decay  constants  varying 
with  distance,  features  that  cannot  be  reconciled  within  the  concept  of 
universality.  Similar  behavior  has  been  found  in  an  earlier  calculation  [9].  The 
rescaled  potential  for  the  Ir-AI  interaction  (dashed  line)  is  qualitatively  similar  to 
the  Al-.AJ  potential  except  that  it  decays  more  rapidly  at  large  separations.  To 
test  the  reliability  of  the  ASED-MO  method,  we  have  calculated  the  AI-Al 
interaction  (I)  using  spin  density  functional  theory  with  both  local  and  nonlocal 
approximations  to  the  exchange  and  correlation  cnerg>'  [10].  As  expected,  we  find 
a  smaller  binding  energy  for  the  Al-Al  system  but  the  rescaled  curve  is  in  very 
good  agreement  with  the  ASED-MO  results.  As  a  final  check  on  the  quality  of 
the  effective  two-body  interaction  (I)  as  determined  by  the  ASED-MO  method. 


we  calculate  the  energy  of  a  large  cluster  of  up  to  211  metal  atoms  by  summing 
nearest  neighbor  interactions  pain-vise  and  minimizing  the  lattice  constant.  For  the 
latter  we  typically  achieve  agreement  wnth  experimental  values  of  within  0.05  A. 

In  tlie  first  set  of  caIculation.s  of  the  interaction  energies  and  the  force 
gradients  between  an  Ir  tip  and  a  planar  fr  and  AI  surface,  we  kept  the  lattice 
spacings  in  both  clusters  at  their  equilibrium  positions,  see  dashed  curves  in 
Fig.  3.  Note  that  the  nominal  gap  width  Z  measures  the  distance  between  the 
apex  atom  and  the  top  atomic  layer  of  the  sample.  The  nominal  gap  width  minus 
a  constant  offset  corresponds  to  the  experimental  tip  excursion  scale  z.  The 
experimental  data  (Fig.  1)  is  reproduced  rather  well  for  the  Ir  surface.  However, 
the  calculated  force  gradient  of  the  Al  surface  deviates  significantly  from  the 
measured  data;  being  in  particular,  one  order  of  magnitude  too  large  (right-hand 
scale  in  Fig.  3b).  Good  qualitative  agreement  is  achieved  by  replacing  the  Ir 
atom  at  the  apex  of  the  tip  by  an  Af  atom  (dotted  curve).  It  is  conceivable  that 
the  tip  accidentally  touched  the  AI  .surface  in  the  course  of  the  experiment  and 
that  some  Al  atoms  were  picked  up  by  the  tip  (note  that  the  tunneling  current 
exhibited  substantial  fluctuations  with  the  AI  samples).  The  measured  force 
gradient  is  slightly  larger  than  that  predicted  by  the  model  calculation.  The 
discrepancy  can  be  explained  by  the  fact  that  the  latter  assumes  a  single  atom  tip 
whereas  such  a  perfect  tip  geometry  was  most  iikely  not  achieved  in  the 
experiment. 

In  the  second  calculation  the  positions  of  the  tip  and  sample  atoms  are 
allowed  to  relax  in  response  to  the  interaction  forces  that  keep  the  base  layers  of 
each  cluster  fixed.  Tip  and  surface  relaxations  along  the  tip  axis  are  shown  in 
Fig.  4  as  solid  and  dashed  curves,  respectively.  Relaxation  causes  the  actual  gap 
width,  defined  as  the  average  separation  between  the  tip  atom  and  the  nearest 
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atoms  in  the  surface,  to  decrease  more  rapidly  than  Z.  Correspondingly  the  force 
gradient  curves  also  change  as  shown  by  the  solid  curves  in  Fig.  3.  Relaxation 
effects  are  particularly  pronounced  for  the  Ir-Ir  system  (Fig.  4a).  At  the  onset  of 
making  contact,  the  apex  atom  on  the  tip  has  moved  towards  the  surface  by  as 
much  as  0.5  A  with  a  smaller  upward  motion  of  the  sunace  atoms  under  the  tip. 


^fn"  uhe  ateni'  "more  bf""fess  fidvehng^'bccweOT'''Cne“  tW^cflisters""tb  '■fbfm'’'a 
metailic  bridge.  This  sudden  change  explains  why  the  tunneling  resistance  drops 
abruptly  at  very  small  gap  widths. 

For  the  Al-coated  Ir  tip  approaching  an  AI  surface,  much  smaller  force 


gradients  are  obtained  as  a  consequence  of  the  - peculiarity  of  the  ,Ai-Al  potential. 


see  Fig.  2.  Correspondingly,  relaxation  of  the  atomic  positions  (Fig.  4b)  is 
substantially  smaller;  it  does  not  exceed  a  tenth  of  an  angstrom.  Note  that  the 
outwards  motion  of  the  apex  and  surface  atoms  is  reversed  when  the  force 
gradient  is  positive.  This  feature  stabilizes  the  gap  to  a  certain  degree.  A  steep 
rise  of  the  atomic  displacements  is  also  obtained.  The  discontinuity  is  not  nearly 
as  dramatic  as  for  the  Ir-lr  system  and  It  occurs  much  closer  to  the  equilibrium 
position  where  atomic  displacements  become  zero  again.  No  evidence  for  a  jump 
to  contact  was  observed  in  the  experiment,  however,  and  this  might  be  explained 
by  the  smallness  of  the  atomic  displacements. 

As  a  final  point  we  should  comment  on  the  fact  that  the  distance  scale  of 
the  theoretical  force  gradients  is  smaller  man  tne  experimental  data  by  aoout  a 

factor  of  1.5.  We  can  see  two  reasons  for  this  discrepancy,  (i)  Our  effective 

nearest  neighbor  interactions  account  well  for  only  the  short-ranged  interactions, 

but  should  be  modified  at  large  separations  by  a  more  slowly  varying  dispersion 

like  force,  (if)  Our  model  calculation  assumes  cIo.se-packed  surfaces  whereas 


polycrystalline  samples  were  investigated  in  the  experiment  and  hence  the  surfaces 
of  the  samples  were  not  perfectly  smooth  on  an  atomic  scale.  Therefore  the  actual 
gap  width  fluctuates  slightly  depending  on  the  lateral  tip  position  which  in  turn 
leads  to  a  broadening  of  the  measured  force  gradient  curves. 

In  summary  qualitative  differences  in  the  adhesion  interaction  of  an  Ir  tip 
with  Ir  and  Ai  surfaces  have  been  observed  experimentally.  They  can  be 
understood  to  arise  (i)  from  small  deviations  from  universality  in  the  effective 
nearest  neighbor  interaction  and  (ii)  from  the  deformation  of  the  lattice  as  the  tip 
approaches.  There  is  a  direct  correlation  between  the  interaction  force  gradient 
and  these  distortions  which  in  turn  detennine  the  formation  of  an  atomic  scale 
contact  Our  model  calculations  provide  a  semi-quantitative  explanation  of  the 
experimental  data.  We  have  discussed  the  reasons  for  small  discrepancies. 
Improvements  can  be  made  by  setting  up  (vastly  .more  complicated  and  less 
transparent)  ab  initio  cluster  calculations  that  account  for  many-body  and 
long-range  interactions. 
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Figure  Captious 


Fig.  1: 

lateraction  force  gradient  VF  and  tunnel  resistance  Rj  versus  tip  displacement  z 
measured  rrith  liie  STM  using  an  Ir  tip  and  (a)  polyci^'stallinc  Ir  (b) 
poiycrystalline  Ai  as  sampie. 


Fig.  2: 
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Fig.  3: 

Force  gradients,  dashed  and  dotted  lines  without  relaxation  and  solid  lines  with 
relaxation,  (a)  Ir  tip  above  an  Ir  surface,  (b)  Ir  tip  with  an  AI  atom  at  its  apex 
above  an  AI  surface;  dotted  line  for  a  clean  Ir  tip  (note  change  in  scale). 

Fig.  4: 

Tip  (solid  lines)  and  snnace  (dashed  lines)  relaxations  for  (a)  the  Ir-Ir  and  (b)  the 
Ir(.A[)-Al  systems. 
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